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Abstract 

We have calculated the evolution of low-mass X-ray binaries that lead to the 
formation of the binary radio pulsars, like PSR J1713+0747. We showed that the 
mass transfer is most likely to be nonconservative, due to unstable disk accretion, 
to account for the mass of PSR J1713+0747, which is close to its initial value. We 
assumed that part of the lost material from the binary may form a circumbinary disk, 
and found that it can significantly influence the mass-transfer processes. We briefly 
discuss the implications of the circumbinary disks on the evolution of low-mass X-ray 
binaries and the formation of low-mass binary pulsars. 
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1. Introduction 

Since the first radio pulsar was discovered by Hewish et al. (1968), the number of pulsars 
has already increased to be more than 1500. In this population only ~ 3% are the members of 
binary systems. Most of the binary radio pulsars are millisecond pulsars (Backer et al. 1982) 
with a He or CO white-dwarf companion. Such binary and millisecond pulsars are thought to 
have been “recycled” from accreting low-mass or intermediate-mass X-ray binaries (L/IMXBs, 
see Bhattacharya, van den Heuvel 1991; Tauris, van den Heuvel 2006 for reviews). Mass 
accretion onto a neutron star induces magnetic held decay, and spins the star up to a short 
period. When mass transfer ceases, the end-point of the evolution is a circular binary containing 
a neutron star visible as a low-held, millisecond radio pulsar, and a white dwarf, the remaining 
He or CO core of the companion. The recent discovery of millisecond accreting pulsars has lent 
strong support to this scenario (van der Klis 2006). 

The mechanisms that drive mass transfer in LMXBs depend on the initial separations 
of the binary components (Bhattacharya, van den Heuvel 1991). In narrow systems with 
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initial orbital periods of P m b < 1 — 2 d, mass transfer is driven by the loss of orbital angular 
momentum via gravitational radiation and/or magnetic braking. Mass transfer in relatively 
wide (P 0 rb ~ 1 — 2 d) LMXBs is driven by the nuclear expansion of the secondary. Systems 
of this kind form a quite homogeneous group whose evolutionary history seems to be well 
understood (Webbink et al. 1983; Taam 1983). Ritter (1999) derived a simple analytical solution 
for the evolution of a close binary with nuclear time-scale driven mass transfer from a giant, 
based on the well-known fact that the luminosity and the radius of a giant scale to a good 
approximation as simple power laws of the mass of the degenerate helium core. It has been 
shown that the orbital binary periods gradually increase during evolution. Once the initial 
orbital period, P or b, is long enough, the mass transfer rate, M, will exceed the Eddington 
accretion rate, Me — 1.5 x lCT 8 M 0 yr _1 , of a neutron star, probably leading to mass loss, and 
the mass transfer becomes nonconservative. 

The 4.57 ms radio pulsar PSR J1713+0747 is in a 67.8 day circular orbit with a low-mass 
white-dwarf companion (Foster et al. 1993). These characteristics indicate that this pulsar 
probably had evolved from a LMXB consisting of a neutron star and a low-mass red giant. 
Recent observations (Splaver et al. 2005) constrain the masses of the pulsar and secondary star 
to be (1.3 ± 0.2) M 0 and (0.28 ± 0.03) M 0 , respectively (68% confidence), implying extensive 
mass loss during the previous LMXB evolution. 

In this paper we explore nonconservative mass transfer in the LMXBs evolution, leading 
to systems like PSR J1713+0747, while taking into account mass loss resulting from unstable 
mass transfer, as well as its feedback to the binary evolution. In section 2 we describe our model 
for the evolution of LMXBs. In section 3 we present the calculated results, and compare them 
with the observational data of PSR J1713+0747 in section 4. We present a brief discussion and 
conclude in section 5. 


2. Model 


In a binary system consisting of a neutron star (of mass Mx) and a low-mass red giant 
(of mass M), the evolutionary expansion of the Roche lobe-filling giant will transfer mass to the 
neutron star through the inner Lagrangian point on the nuclear time scale. The mass transfer 
rate can be described by the following equation (e.g., Ritter 1999): 


Mi 


M 


I , * > (!) 

M 

where £ e = (<9 In R/d In M) eq is the thermal equilibrium mass radius exponent and = 
dlnR^/dlnM is the mass radius exponent of the Roche lobe radius of the giant (Soberman et 
al. 1997); R and are the radii of the giant and of the Roche lobe, respectively. 
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2.1. Nuclear Evolution of the Giant 


For a lower giant-branch star, the stellar luminosity, L, and radius, R, are uniquely 
determined by the mass of its degenerate helium core, M c , independent of the hydrogen envelope 


mass (Refsdal, Weigert 1970), 


L = L e exp (Sep?/*), 

(2) 

R = P 0 exp(Xqy*), 

(3) 


where y = ln(4 M c /M & ) and i — 0, 1, 2, 3. The parameters a, and c* in the above equations for 
the solar chemical composition can be found in Webbink (1975). 

From equation (3), the change in the radius of the donor star is related to the growth 
rate, M c , of the helium core by 

/ R\ 9n M c 

w„ =(ci 2S 3B) ^' (4) 

Because the stellar luminosity is directly related to the growth in the core mass due to 
shell hydrogen burning, 

L = Xe h M c , (5) 


where X is the initial hydrogen mass fraction in the shell, and £u ~ 6 x 10 18 erg g -1 is the 

released energy rate by the CNO cycle (Webbink et al. 1983). 

Since the giant remains in thermal equilibrium during the mass transfer, and its radius is 
only determined by the core mass, independent of total mass, M, csO (Soberman et al. 1997). 
We then have 


M 


ci + 2 c 2 y + 3 c 3 y 2 
£r 



)M C . 


( 6 ) 


2.2. Mass and Angular Momentum Loss 


The orbital evolution is determined by the angular-momentum loss rate for a binary 
system possessing mass communion. The total angular momentum of a system with a circular 
orbit is 

T _ 2 MM X 27t 
' =f ' M + MxP 0 r b 

where a = [G{M + Mx)Po rb /(47r 2 )] 1 / 3 is the binary separation; G and P or b are the gravitational 
constant and the orbital period of system, respectively. We neglect the spin angular momentum 
of the components because it is small compared the total orbital angular momentum of the 
system, and consider the following processes for the mass and orbital angular momentum loss 
from the binary system. 

1. isotropic wind 

We consider the mass and angular momentum loss during nonconservative mass transfer 
in the evolution of LMXBs. If the magnitude of the secular mass transfer rate, —M, is greater 
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than the Eddington accretion rate, Me, we assume that the neutron star accretes a fraction 
(1 — a) of the transferred mass, and the remaining fraction, ct, is ejected in the vicinity of the 
neutron star (a = 0 if |M| < Me, else a = 1 + Me/M). 

Mass loss can also occur when the mass transfer is sufficiently low. It is well known that 
disk accretion in LMXBs is thermally and viscously unstable if the mass-transfer rate is less 
than a critical value, M cr (van Paradijs 1996; Dubus et al. 1999) 


M cr ~ 3.2 x 10 9 M q yr 
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When an accretion-disk instability occurs, the accreting neutron star will become a transient 
X-ray source, experiencing outbursts separated by long quiescent intervals. For the accretion 
behavior during outbursts, we adopt the following prescription suggested by Portegies Zwart 
et al. (2004). The accretion rate, M x , first reaches a peak value of 2Me, and then decays with 
an exponential timescale of td = 6d x min(l, P or b/10hr). The relation of the recurrence time, 
t r , and the total energy, E, in an outburst satisfies t T = E/(— O.lMc 2 ), where log(P/erg) = 
45 + log(P orb /d). Since the accretion rate of a neutron star is M x = 2M E exp(— t/t d ), we can 
calculate the accreted mass of the neutron star within the recurrence time, t T , to be 


A M x = 2M E l d [l -exp(-« r /( d )]. 


(9) 


Accordingly, the fraction a = 1 +AM x /(Mt r ) of the transferred mass is assumed to be lost 
from the system. 

The material probably leaves the binary system in the form of winds, outflows, or jets. 
We assume that a fraction of (1 — 5) of the lost matter is ejected in the vicinity of the neutron 
star, carrying away a specific orbital angular momentum of the neutron star, 


Jw M t M x ’ 

where Mp is the total mass of the system. 
2. CB disk 


( 10 ) 


Spruit and Taarn (2001) suggested that a Keplerian, circumbinary (CB) disk could be 
formed as the result of mass outflow from an accreting star or accretion disk. Tidal torques are 
put on the CB disk, and carry away angular momentum from the binary orbiting inside it. We 
assume that the other part (with a fraction of h) of mass loss forms a CB disk, extracting the 
orbital angular momentum from the binary system. 

At the inner edge, r ; , of the disk, the viscous torque exerted by the CB disk on the 
binary can be shown to be (Spruit, Taarn 2001; Taam, Spruit 2001) 

r = 7(jjA«(/-) 1/2 , (ii) 

-‘orb ^vi 

where, y 2 = r-Ja , a is the binary separation, and t vi the viscous timescale at the inner edge of 
the CB disk, 
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where V\ is the viscosity at the inner edge of the CB disk. We estimate zq using the standard a 
prescription (Shakura, Sunyaev 1973), 

R = a ss c sHi, (13) 


where c s and 77; are the sound speed and the height of the disk at the inner edge of disk, 
respectively. In the following calculations we set the viscosity parameter, ass — 0.001, and 
assume that the disk was hydrostatically supported and geometrically thin with H-Jr^ ~ 0.03 
(Belle et al. 2004). The sound speed can be obtained from the equation of vertical hydrostatic 
equilibrium, 


c s ~ fl K 77j, (14) 

where fix = (GMx/rf) 1 / 2 is the Keplerian angular velocity at r\. 

3. the gravitational and baryonic mass of the neutron star 

For neutron stars, one should consider the discrepancy between its baryonic mass, M b , 
and the gravitational mass, M g . By differentiating the binding energy equation (Lattimer, 
Yahil 1989), 

M b = M g + O.O84M 0 (^) 2 , (15) 

M 0 

with respect to t, the relation between the growing rate of the gravitational and of the baryonic 
masses is 


Mx = M g = (l-/?)M b , (16) 

where (3 — 1 — (1 + 0.168Mx/Mq)^ 1 . We assume that the fraction f3 of the accreted baryonic 
mass disappears along with the specific angular momentum of the neutron star. 

2.3. Computing 


Considering various kinds of mass loss discussed in the above subsection, we can write 
the rate of change of the orbital angular momentum of the binary to be 


- = [(1 - 5 ) ——— 
J J M t M x 


+ 77—laM + 

y m t m x 


M 


(1 -a)PM, 


(17) 


where r) = S(t/t vi yG is the efficiency factor of transfer mass removing angular momentum from 
the system through the CB disk; p = MMx/Mp is the reduced mass of the binary. The first, 
second, and third terms on the right-hand side of equation (17) represent the change in the 
orbital angular momentum due to mass loss through ejected outflows, the CB disk and the 
baryonic mass loss, respectively. 

For M < 0.8Mx, the Roche lobe radius of the donor can be approximately denoted by 
(Paczynski 1971) 
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R l = 0.462a —C- 

\l + q 

where q = M/M x . We then obtain £ R combining equations (7), (17), and (18): 

£r — — g + 2(1 — a)(l — {3)q + -(a + {3 — a/3) ^ + 

2[a(l - 5) +/3(1 - a)] + 2 7?/«(l + ?)■ (19) 

When a = /3 = 0, equation (19) recovers to the standard conservative mass-transfer case, 

£ R = -| + 2g. (20) 

When {3 = 5 = 0, equation (19) becomes 

& = -5 + 2(1 -«)«+(! + 29 )^ ( 2 !) 

for nonconservative mass transfer with only outflows indued (Li, Wang 1998). 

3. Numerical Results 

We adopted the semi-analytical method by Webbink et al. (1983) to calculate the 
evolutionary sequences of LMXBs. We set the initial mass of the donor M = 1.0M 0 (with solar 
chemical composition X=0.7, Y=0.28) and the core mass in the range of O.15M 0 <M C <0.35M 0 , 
the initial mass of the accreting neutron star being M x = 1.4M 0 . 

For stable mass transfer, we let M x = — (1 — a)(l — {3)M when — M > M R . When the 
accretion disk becomes unstable, we calculate AM x according to equation (8). For the radius r\ 
of the inner edge and the mass feeding rata 5 of CB disk, we took r-Ja = y 2 = 1.7 (Artymowicz, 
Lubow 1994) and 5 = 0, 0.01 and 0.02. We stopped the calculations when the mass of the 
hydrogen envelope of the giant was reduced to be ~ 2% of the core mass. 

Figure 1 shows the calculated £ R against the mass ratio, q = M/M x , for the initial 
M c = O.25M 0 and 5 = 0, 0.01, 0.02. One can see that the larger is 5, the smaller is — £ R , and 
hence the higher the mass transfer rates. This can also be seen in Figure 2, which compares 
the detailed evolutions of the the giant mass, the orbital period, the neutron star mass, and 
the mass-transfer rate for different values of 5. Obviously, the existence of the CB disk can 
accelerate the evolutionary processes quite efficiently. If 5 is further increased, the orbital period 
will decrease, finally leading to runaway mass transfer. The variation of M is also much larger 
if the CB disk effect is included. Because of unstable mass transfer induced by the accretion- 
disk instability, the efficiency of neutron-star accretion is quite low. The resulting averaged 
accretion rates of the neutron star, M x = AM x /t r , are shown in Figure 3. 
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4. Application to PSR J1713+0747 


We have performed many calculations for LMXBs evolution with different initial pa¬ 
rameters in order to fit the observed data (P orb , M x , and M) of PSR J1713+0747. We found 
reasonable results when the initial core mass was O.258M 0 , O.246M 0 , O.2O5M 0 for 8 = 0.02, 
0 .01, 0, respectively, while conservative evolution failed to produce results compatible with all 
of the measured data simultaneously. Table 1 summarizes the results of our calculations. It 
seems that both models with and without a CB disk can well fit the observed data for this 
source, if unstable mass transfer is taken into account. However, there exist considerable dif¬ 
ferences in the initial binary parameters and the evolutionary timescales, suggesting that the 
CB disks may have important implications for the LMXB evolution (see below). 

Previous investigations suggested a simple relation between the orbital period, P orb , and 
the mass, M, of the white dwarf for low-mass binary radio pulsars. Assuming the mass-feeding 
rate to be 5 = 0.002 (a larger 5 possibly cause infinite mass transfer rate when initial core mass 
M c = O.15M 0 , see Figure 1) of the CB disk, we calculated 21 LMXBs evolution sequences, 
where the initial core mass of giant ranged from 0.15 to O.35M 0 . The relation between P orb 
and M when evolution ceases for O.15M 0 < M ci < O.35M 0 was plotted with the solid curve 
in Figure 4; the dot and dashed curves correspond to the relations obtained by Rappaport et 
al. (1995) and Tauris and Savonije (1999). It can be seen that there is a power-law relation 
between P orb and M with P orb (d) oc (M/M 0 ) 6 - 55 . 

5. Discussion and Summary 

The mass measured for PSR J1713+0747 indicates that the neutron star had accreted 
< 0.1 Mq mass during its previous LMXB evolution, if formed with an initial mass of 1.4 M 0 . 
We suggest that unstable mass transfer in the accretion disk may lead to extensive mass loss in 
the evolution of relatively wide LMXBs. The mass accreted by the neutron star is considerably 
smaller than in the stable accretion case (see also Li, Wang 1998). 

The evolutionary path of the binary also depends on the form of mass and angular 
momentum loss when mass transfer is nonconservative. Traditionally the excess material is 
assumed to leave the binary with a specific orbital angular momentum of either the neutron 
star or the binary. However, based on arguments proposed by Spruit and Taam (2001) and 
Taarn and Spruit (2001), we considered the case with a CB disk originating from part of the mass 
outflow. Actually, radio observations of super-Eddington accretion systems SS 433 (Blundell et 
al. 2001) and Cygnux X-3 (Miller-Jones et al. 2004) have shown evidence of equatorial, disk-like 
outflows. 

Our calculations show that the CB disk can remove the orbital angular momentum of 
the binary at a rate significantly higher than simple mass outflows. This may have important 
implications on the LMXBs evolution and its relation with low-mass binary pulsars. Recent 
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binary population synthesis studies (Podsiadlowski et al. 2002; Pfalil et al. 2003) fail to 
produce enough luminous LMXBs, as observed. This problem seems to be related to the 
discrepancy between the birthrates of LMXBs and low-mass binary pulsars (e.g. Kulkarni, 
Narayan 1988), the solution of which may lie in an unknown mechanism for angular-momentum 
loss to speed up the mass transfer. Moreover, it has been pointed out by Sills et al. (2000) 
that the traditional magnetic braking, which has served as the basis for the evolution of 
LMXBs and cataclysmic variables (CVs), is significantly reduced. The CB disk may provide an 
efficient mechanism to drain the orbital angular momentum, accelerating the binary evolution 
and enhancing the mass-transfer rates. Dubus et al. (2002) found that the spectral energy 
distributions expected of CB disks in CVs can dominate the emission from the donor star 
and the accretion disk of the white dwarf at wavelengths £ 3 n m. At longer wavelengths 
the relative contribution from the CB disk to the total emission from the system increases. 
Since LMXBs and CVs have similar binary parameters, searching for (mid- and far-) infrared 
emission from the CB disks will be of great importance for the evolution of both LMXBs and 
CVs. 
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References 

Artymowicz, P., & Lubow, S. H. 1994, ApJ, 421, 651 

Backer, D. C., Kulkarni, S. R., Heiles, C., Davis, M. M., & Goss, W. M. 1982, Nature, 300, 615 
Belle, K. E., Sanghi, N., Howell, S. B., Holberg, J. B., &; Williams, P. T. 2004, AJ, 128,448 
Bhattacharya, D., & van den Heuvel, E. P. J. 1991, Phys. Rep., 203, 1 

Blundell, K. M., Mioduszewski, A. J., Muxlow, T. W. B., Podsiadlowski, P., & Rupen, M. P. 2001, 
ApJ, 562, L79 

Dubus, G., Lasota, J. -P., Hameury, J. -M., & Charles, P. 1999, MNRAS, 303, 139 
Dubus, G., Taam, R. E., & Spruit, H. C. 2002, ApJ, 569, 395 
Foster, R. S., Wolszczan, A., & Carnilo, F. 1993, ApJ, 410, L91 

Hewish, A., Bell, S. J., Pilkington, J. D., Scott, P. F., & Collins, R. A. 1968, Nature, 217, 709 
Kulkarni, S., & Narayan, R. 1988, ApJ, 335, 755 
Lattimer, J. M., & Yahil, A. 1989, ApJ, 340, 426 
Li, X.-D., & Wang, Z.-R. 1998, ApJ, 500, 935 

Miller-Jones, James C. A., Blundell, K. M., Rupen, M. P., Mioduszewski, A. J., Duffy, P., & Beasley, 
A. J. 2004, ApJ, 600, 368 
Paczyriski, B. 1971, ARA&A, 9,183 

Pfahl, E., Rappaport, S., & Podsiadlowski, Ph. 2003, ApJ, 597, 1036 
Podsiadlowski, Ph., Rappaport, S., & Pfahl, E. D. 2002, ApJ, 565, 1107 
Portegies Zwart, S. F., Dewi, J. &; Maccarone, T. 2004, MNRAS, 355, 413 


Rappaport, S., Podsiadlowski, Ph., Joss, P. C., Distefano, R., & Han, Z. 1995, MNRAS, 273, 731 

Refsdal, S., & Weigert, A. 1970, A&A, 6, 426 

Ritter, H. 1999, MNRAS, 309, 360 

Shakura, N. I., & Sunyaev, R. A. 1973, A&A, 24, 337 

Sills, A., Pinsonneault, M. H., & Terndrup, D. M. 2000, ApJ, 540, 489 

Soberman, G. E., Phinney, E.S., <fc van den Heuvel, E. P. J. 1997, A&A, 327, 620 

Splaver, E. M., Nice, D. J., Stairs, I. H., Lommen, A. N., & Backer, D. C. 2005, ApJ, 620, 405 

Spruit, H. C., & Taam, R. E. 2001, ApJ, 548, 900 

Taam, R. E. 1983, ApJ, 270, 694 

Taam, R. E., & Spruit, H. C. 2001, ApJ, 561, 329 

Tauris, T. M., & Savonije, G. J. 1999, A&A, 350, 928 

Tauris, T.M., & van den Heuvel, E.P.J. 2006, in Compact Stellar X-Ray Sources, ed. W. H. G. Lewin 
&; M. van der Klis (Cambridge: Cambridge Univ. Press) in press (astro-ph/0303456) 
van der Klis, M. 2006, to appear in Compact stellar X-ray sources, ed. W. H. G. Lewin & M. van der 
Klis (Cambridge: Cambridge University Press) in press (astro-ph/0410551) 
van Paradijs, J. 1996, ApJ, 464, L139 
Webbink, R. F. 1975, MNRAS, 171, 555 

Webbink, R. F., Rappaport, S., Sz Savonije, G. J. 1983, ApJ, 270, 678 


9 


Table 1 . Calculated results for different evolutionary models. * 


Model 

(<5,a,/3) 

Me (to) 
(Mq) 

tf — to 

(Myr) 

log M(tf) 
{Mq yr" 1 ) 

R(tf) 

(Re) 

-Porb(tf) 

(days) 

Mc(t f ) 

(Mq) 

Mx(tf) 

(Mq) 

,5 = 0.02 

0.258 

23.9 

-7.95 

21.38 

68.36 

0.284 

1.402 

5 = 0.01 

0.246 

47.1 

-8.15 

21.20 

67.71 

0.283 

1.403 

O 

II 

0.205 

267.5 

-8.406 

21.67 

69.78 

0.284 

1.408 

s 

II 

"Co 

II 

o 

0.215 

179.00 

-8.400 

20.99 

66.70 

0.283 

2.117 


* to and tf denote the ages at the beginning and end of the mass transfer, respectively. 



Fig. 1. Plot of vs. q for LMXBs with an initial core mass of M c = 0.25Mq. The dot, dashed, and 
solid curves correspond to <5 = 0, 0.01, and 0.02, respectively. 
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Fig. 2. Evolutions of an LMXB with an initial donor mass of M = 1.0M 0 and a core mass of M c = 0.25 Mq. 
The dot, dashed, and solid curves correspond to 5 = 0, 0.01, and 0.02, respectively. 
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Fig. 3. Evolution of the accretion rate of the neutron star (mean accretion rate in a recurrence time f r ). 
The dot, dashed, and solid curves correspond to 5 = 0, 0.01, and 0.02, respectively. 
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Fig. 4. Predicted relation between P or b and M for low-mass binary radio pulsars. The dashed, dot, and 
solid curves correspond to the relations obtained by Tauris and Savonije (1999), Rappaport et al. (1995), 
and this work, respectively. 
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